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CD I The local and global thermal phase structure for asymptotically anti-de Sitter black 

holes charged under an abelian gauge group, with both Gauss-Bonnet and quartic field 
strength corrections, is mapped out for all parameter space. We work in the grand 
canonical ensemble where the external electric potential is held fixed. The analysis 
^ ' is performed in an arbitrary number of dimensions, for all three possible horizon 

■ topologies - spherical, fiat or hyperbolic. For spherical horizons, new metastable 
configurations are exhibited both for the pure Gauss-Bonnet theory as well as the 
pure higher derivative gauge theory and combinations thereof. In the pure Gauss- 
Bonnet theory with negative coefficient and five or more spatial dimensions, two 
locally thermally stable black hole solutions are found for a given temperature. Either 
one or both of them may be thermally favored over the anti-de Sitter vacuum - 
corresponding to a single or a double decay channel for the metastable black hole. 
Similar metastable configurations are uncovered for the theory with pure quartic field 
strength corrections, as well combinations of the two types of corrections, in three 
or more spatial dimensions. Finally, a secondary Hawking-Page transition between 
the smaller thermally favored black hole and thermal anti-de Sitter space is observed 
when both corrections are turned on and their couplings are both positive. 
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1 Introduction 

A common feature of UV completions of our current theoretical framework, the Ein- 
stein action coupled to the matter content of the Standard Model, are higher deriva- 
tive corrections to both the gravitational and matter low energy effective actions 
[B 121 El HI E] . These corrections are suppressed by the scale where the new effective 
description becomes relevant, which may typically be the SUSY breaking scale or 
the Planck scale to give but two examples. It is therefore important to study such 
higher derivative corrections in order to better understand the properties of the low 
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energy effective actions of possible UV completions to the Standard Model coupled 
to Einstein gravity. 

Of the large variety of possible gravitational corrections, there arises a particularly 
interesting case, where the additional term is given by the Gauss-Bonnet correction 
[B m E]. In four dimensions, this is a topologically invariant term that is analogous 
to the Euler number in two dimensions. Its variation is a total derivative and thus 
does not affect the classical bulk equations of motion. Upon studying gravitational 
perturbations of theories modified by higher derivative terms in any dimension, it 
has been discovered that only the Gauss-Bonnet term does not give rise to ghost-like 
propagating modes with modified kinetic terms [H [2] at the level of Ricci squared 
corrections. The Gauss-Bonnet term, for example, is the a' order correction 
to the Einstein action in the ten dimensional heterotic string theory [3 [H [9], and 
may arise in the lower dimensional actions of various compactified string theorie^l] 
[10 \ i m [12]. Furthermore, it can also result from a' corrections to the bosonic sector 
of the M = 2 gauged supergravity that arises in the case of spinning D3-branes in type 
IIB superstring theory or spinning M2-branes in M-theory. In the dual conformal field 
theory, these corrections would be of order 1/\/A in the strong coupling expansion. 

Another endless source of fascination has been the discovery that black holes 
have a very natural thermodynamic description [13^ I14j . In particular, they emit a 
thermal spectrum and have an entropy related to their surface area. Of significant 
interest are black holes in asymptotically anti-de Sitter space, since these black holes 
have the remarkable feature of being thermodynamically favored for high enough 
temperatures, as was shown in the classic paper |15j . With the advent of the anti-de 
Sitter/conformal field theory correspondence |161 [T71 118j . such anti-de Sitter black 
holes have become even more interesting, since their thermodynamic structure is 
related to that of the dual conformal theory [l9l[20l[2ll[22l[23l[2ll[25l[26l[27]. 
Specifically, the high temperature black hole phase is dual to a deconfining phase, 
whereas the low temperature thermal anti-de Sitter phase is dual to a confining phase 
in the boundary theory [25] . 

In this paper we will study asymptotically anti-de Sitter black hole solutions, 
which are electrically charged under a gauged U{\). Our gravitational action is 
corrected by a Gauss-Bonnet term, and we further include an correction to the 
electromagnetic field strength in the matter action. These corrections are a subset of 
all possible a' corrections to the bosonic sector of = 2 supergravity in four or five 

^In pure IIB string theory on AdSs x , there are no corrections [3] ID [5], but they may arise 
when we compactify on less symmetric manifolds. 
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spacetime dimensions [12] , and also a subset of a' corrections for higher dimensional 
anti-de Sitter space with a gauged U{1). We obtain explicit black hole solutions and 
study their respective global and local thermodynamic stability in the grand canonical 
ensemble, where the external electric potential is held fixed. The thermodynamic 
phase structure we obtain should, in principle, be related to the thermodynamic phase 
structure of the dual boundary gauge theory, with a constant chemical potential. We 
will build our results from the well studied cases of Reissner-Nordstrom anti-de Sitter 
black holes [IHl El [Ml [29] and Gauss-Bonnet black holes [261 [13 [SQl [3ll [32l [33] , and 
generalize to the Gauss-Bonnet-F^ black hole. 

There are various features that are novel to this discussion. As far as we know, 
the thermodynamics and thermal phase structure of the Gauss-Bonnet-F^ black hole 
has never been studied before. It is thus of interest to study the thermodynamics of 
three types of limits in order to understand what the features of the thermodynamic 
structure are attributed to. The first system we will explore is the charged Gauss- 
Bonnet black hole with no F'^ corrections. Then we will study the -F^ black hole with 
no Gauss-Bonnet corrections and finally we will obtain the thermodynamic structure 
of the black hole with both corrections turned on. Furthermore, we will study black 
holes with hyperbolic, spherical and flat spatial horizons [Ml [SS [361 [23 [Ml [39] - 
known as topological black holes - and express our results for an arbitrary number of 
dimensions. 

Due to the various systems we are studying, we will organize the paper in three 
parts. The first part of the paper will discuss the action of our theory and the explicit 
solutions to the ansatz. We will further identify the parameter space giving rise to 
actual black holes. In the second part, we will review the Gauss-Bonnet black hole 
with no F'^ corrections. The results will be presented in a parametrization that is 
most convenient when studying the F^ black hole for the grand canonical ensemble, 
which is the ensemble most relevant to holography. In the third part of the paper 
we present the thermodynamic structure of the Gauss-Bonnet- i^^ black hole, again 
for the grand canonical ensemble. We conclude the paper with a brief discussion and 
summary of our results. 
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PART I - The Black Hole Solution 



2 Framework 

We will now develop the basic setup of the gravitational and matter actions. This 
will be the goal of the section below, where we derive the equations of motion given a 
stationary spherically symmetric ansatz for the metric and the assumption of purely 
electrically charged matter. 

2.1 i?^ Corrections and the Gauss-Bonnet Term 

Our story begins with the action for Einstein-Maxwell gravity with negative cosmo- 
logical constant and higher derivative corrections. In our conventions we have the 
following action: 

Itot — Igrav ~l~ Imatten (^•-'^) 

where 

Igrav = d'^+'x^/^ (^^R - A + aR^ + bR'^-'R^, + cR'^-'P'^R^.p^^ . (2.2) 

As we shall see, the Hawking-Gibbons boundary term will not play a role in our 
discussion and so we do not include it in the action. 

From the point of view of a stringy UV completion these R^ corrections would 
have coefficients of order a' . A particular combination of the coefficients, given by 
a = c = —b/4:, gives rise to the Gauss-Bonnet term. In four dimensions, this term 
gives rise to a total derivative upon variation of the action with respect to the metric 
and thus is a topological invariant for four dimensions. It is interesting to note that 
of all ten-dimensional string theories, it is only the heterotic string theory that picks 
up R^ corrections, and they always come in the form of the Gauss-Bonnet invariant. 
Furthermore, it is only the Gauss-Bonnet invariant that does not give rise to ghost- 
like propagation for gravitational perturbations [2j. For such reasons we will work 
exclusively with the Gauss-Bonnet correction throughout the rest of the discussion. 

Having specified our gravitational action, we can obtain the equations of motion 
for some arbitrary matter content. We will use the following convention for our 
stress-energy tensor: 

1 ^Imatter 

Upon variation of the action Igrav with respect to g^i, we obtain the equations of 



rp ^ '-'-^■matter r,\ 
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motion for our gravitational theory with Gauss-Bonnet corrections: 



-2 i? - A + C _ ARf.^R"" + Rf^upaR"""") 



H — 2^1^'^ ~^ ^ ('^RRfiiy — ^RiipRy — ^RppvcjR^^'^ + 2R^^'^ R^ 



^uXpcr 



(2.4) 



For the benefit of the reader we have included a rather lengthy derivation of this and 
subsequent results, which demand extensive computations, in appendix [XI In what 
follows we will search exclusively for spherically symmetric solutions. In particular, 
we will be using the stationary spherically symmetric ansatz 



d-l 



»J=1 



(2.5) 



where we have that gij is an Einstein metric whose Ricci tensor satisfies Rij = kgij, 
and furthermore, we rescale the such that det[5jj] = 1. Many of the steps required 
to obtain the equations of motion are given in appendix [Bl We reproduce here only 
the tt-equation and rr-equation which will be of more immediate use in the discussion 
that follows. The tt-equation of motion is given by: 



(d-l) 



2r^ 



+ e 



(d-l) (d- 3)c 



1 (d-4)^ 

2 (d-2)r4^ 



2r2 



1 



-2u 



A 



2kX' {d-4)k 



(2.6) 



and the rr-equation of motion is given by: 



d-l 



-2X 



{d- l)(d-3)c 



k 

1 _^^{d-A)k^ 



+ le-A 



+ e 



-4A 



2 (d - 2) r4 

+ {d-2) e 



2ku' (d-4)/c 

.^^f2u' l (d-4) 
.3^2 r4 



(2.7) 



Our next task will be to obtain the equations of motion and stress-energy tensor 
following from the matter Lagrangian. 
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2.2 Field Strength Corrections and the Term 

Our matter content is described by the usual Einstein-Maxwell matter Lagrangian 
with a quartic field strength correction, given as follows: 



1 

47 



i matter — A 2 i ^ ~QF p^yF^ 



+ y d'^^^x^g [ci {FpyFn^ + ^^Fp.F-^F^.FP^ 



(2i 



Once again from the stringy point of view, the F'^ term would be of order a', which is 
the same order as the Gauss-Bonnet term. In fact, it is more general to include other 
a' terms of the forms F^R and (V-F)^ (see for example \yi\ I40j). Furthermore, there 
is a potential Chern-Simons term of the form C5(j4^) = e^^f"'^ Ap,FypFfjT which could 
be included. For purely electric solutions, however, this term will not contribute to our 
equations of motiorJ§. In any case, it is not our purpose here to obtain the most general 
OL corrections to the Maxwell-Einstein theory, but to study the thermodynamics of 
a theory with both Gauss-Bonnet and gauge field strength corrections in asymptotic 
anti-de Sitter spacfl At the linearized level, however, it was shown in [H] that field 
redefinitions of the metric and gauge field allow us to eliminate four of the eight 
possible 0(0') terms in the effective Lagrangian, at least for d = 4 spatial dimensions. 

Having established our matter Lagrangian we can compute the stress energy- 
tensor by applying on Imatter- 



+ ci { ^^g'^" [FP^Fp^f - AFxPF^^FP'^Fp^ 



^More precisely, we will be working with only Ftr{r) non-zero from which it is evident that there 
is no non-vanishing Chern-Simons term. 

^In fact as pointed out in [40], we can take an example of this situation to be the ten-dimensional 
Lagrangian of the -Eg x Eg or 50(32) heterotic string theory. The Eg x Eg or 50(32) gauge group 
has a U{1) subgroup with field strength -F^j^. Turning off all other gauge fields and higher forms and 
keeping the dilaton 00 constant leaves us with the ten dimensional effective quartic order Lagrangian 
[ZllSlin] (albeit lacking a cosmological constant): 



R p p*^" J 



r — R ^ (? p*^" _L I n ntf^p" A T3 pA"^ j„ r?2\ 



64 

where h = e""^"'^. 
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+ C2 (^^g^'-F^PFp^F^^F^^ - 4F^^Fp,F'^-F/) . (2.9) 

We proceed to obtain the equations of motion for the Maxweh field strength. These 
are a modified form of the Maxwell equations in a curved background: 



g 



-g ( -^F'"' + 8ci {FP^'Fp^) F^'^ + ScsF^^'F^^F'^'^ 



0. (2.10) 



We note that the equations will depend only on the combination (2ci + C2), which 
we denote by e from now on: 

e = 2ci + C2. (2.11) 

Our next task is to solve these equations for the spherically symmetric case, and we 
will further assume that only Ftr is non-vanishing. For convenience we define the 
following: 

Ftr{r)=e''+^f{r), (2.12) 
such that the non-vanishing components of the stress-energy tensor become: 

Ag 



T« = e-2-(^/(r)^ + 3e/(r)^), (2.13) 



= -e-2A(^-L/(,)2 + 3e/(r)^j, (2.14) 

= ^(4^/WVe/(r)^). (2.15) 

Inserting the ansatz (|2.12|) into the modified Maxwell equations gives rise to the 
following equation for /(r): 

8sf{rf + ^f{r)--^ = 0. (2.16) 

This cubic equation can be solved exactly, and will have either three or one real 
solutions, from which we have to choose the physical one. By physical we simply 
mean that solution which leads to the usual Coulomb's law when e approaches zero 
and which does not diverge as r approaches infinity. The general solution to the cubic 
is given as follows: 

27rn ^ ^ ^ \ 
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where n G {0, 1, 2} and we have defined: 



A(r) 



1/3 



(2.18) 



When e > the physical solution corresponds to the only real solution which is given 
by n = 0. When e < the physical solution is given by n = 1. This solution is 
real only for radii larger than a critical value, rmin, obtained by requiring that the 
expression in the root within A(r) is positive. In particular, we find that 



(2.19) 



For both of the physical solutions there are two complementary expansions to con- 
sider. For X = {8g^Q'^e)/r'^^'^^^^ ^ 1 we can expand as follows: 



fir) 



(2.20) 



We immediately recognize the first term as Coulomb's law which is subsequently 




f(e=0) 



/„(e=10-^) 



/,(e=-10-') 




(a) 



(b) 



Fig. 1: (a) The physical solution for e = 0, e > and e < 0, and (b) the validity of 
the X 1 and X ^ 1 expansions for /(r). 



corrected by e effects. When x ^ = '^^^'^ /i^g^Q'^s) <C 1 then the above expansion 
breaks down and we should consider the following expansion: 



fir) 



Q 



8er" 



1/3 



X 



-1/3 



+ o{ix 



(2.21) 
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The physical solutions and the validity of the two expansions is displayed in figures 
[T]^a) and (b). When e < 0, the region < 1 is excluded due to the physical bound 
(j2.19p . so we should always consider the x expansion in this case. In what follows x 
is always evaluated at r = rn, the horizon radius. Finally, when Q = the physical 
solution corresponds to /(r) = 0, since the other two blow up as e tends to zero. 

Having established our basic conventions and formalism for the higher derivative 
corrections to the Einstein-Maxwell theory, we are ready to find the corresponding 
black hole solutions. This will be the main goal of the next section. 



3 The Gauss-Bonnet-F^ Topological AdS^+i Black Hole 

We have obtained the dynamical equations for our action with a given matter content 
and we are now in the position to obtain an explicit solution for our ansatz. We will 
identify this solution as a static charged black hole in asymptotic anti-de Sitter space, 
and we evaluate its Hawking temperature and horizon radius. 

3.1 Exact Solution 

In order to obtain the black hole solutions to the action we have specified, we can 
take advantage of some of the equations of motion we derived in the previous section. 
In particular, we find the following useful relation: 



= e^'^T" + e^^T"' 



(3.1) 



The following result is immediately yielded: 



ly' + X' = 0. 



(3.2) 



This allows us to set v = —A up to a constant shift in i^, which can always be 
reabsorbed into our definition of time. Having obtained this result, some algebra on 
the tt-equation of motion given in (j2.6p gives rise to the equation: 



(d-1) 



d_ 

dr 



+ {d-l){d-3)c 



{d - 4) k^r 



2^ ( e-2AA;r^"4 



dr 



-^(e-"(.-2)/- 
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+ 2r' 



It is a simple task to integrate over r and find: 



d-i 



{rf + 2>ef (r) 



0. (3.3) 



id -I) 



2Xd-2 



r 



+ 



1 



{d- l)(d-3)c 



d-2 



-kr 



d-2 



(d-2) 



{d - 2) r' 



d-A 



where we have defined: 



/(r) 



drr 



d-l 



—firy + 3ef{r) 



+ /(r) = /i, (3.4) 



(3.5) 



and /i is an integration constant, which we will soon relate to the mass of the black 
hole. We select the integration constant so that the limit of I{r) at infinity vanishes, in 
order to make connection with the usual Reissner-Nordstrom solution in the absence 
of corrections. The above integral /(r) will play a crucial role for the remainder of 
the paper, for it encompasses the e dependence of the metric in a highly non-trivial 
way. Its expansion for x ^ ^ will be of use in later sections and is given by: 



lirn) 



2{d-2)r%~ 



d-2 



(d-2) 6(d-2) 2 



(3.6) 



Notice that (j3.4p is a quadratic equation in e which can be solved quite simply 
to give our metric component: 



-2A 



d-2 



2k'^D 



k ^{d-l)r^ _^ K'^r^ ^{I{r)- ii) 



4Z)2 



jjrd-4 ' 



(3.7) 



where 



D = c{d-l){d-2){d-3) 
2 _ (d-lf ADA 



(3.8) 
(3.9) 



For particular regions of parameter space, the solution we have obtained contains a 
physical singularity at the origin hidden by event horizons and is thus a black hole. 
We proceed to examine the relevant parameter space. 
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c> 



Let us discuss first the case c > for our solution (j3.7p . We take A; > for the 
discussion in this subsection, and leave the case of k < with c > for the next 
subsection. Of the two roots in (|3.7|) only the negative one reduces to the usual 
anti-de Sitter Reissner-Nordstrom solution as c and e go to zero. This will be the one 
we will work with from now on for c > 0. It has actually been shown for the e = 
case that the positive root solution (with ck > 0) has gravitons that have ghost-like 
propagation [U and we take this as further reason to stick to the negative root. 

It is also clear from the form of the metric that arbitrarily large values of c > 
are not acceptable. In fact, it is clearly seen that as r — > oo the values of c that give 
rise to a real metric must satisfy the following bound: 



If c is above this bound our metric becomes complex for the particular ansatz we 
have chosen and is thus rendered unphysical. It was, however, pointed out in [31] 
that exactly where the asymptotically anti-de Sitter metric breaks down the de Sitter 
metric with positive A yields an acceptable solution. We will not consider the im- 
plications of this observation and simply end our parameter space where the bound 
(j3.10p is saturated. 



For the case of c < with k > our solution (j3.7p tends to the usual anti-de Sitter 
Reissner-Nordstrom black hole as c — > for the positive root. Furthermore, we must 
be particularly careful as to when the positive root solution of (j3.7p gives rise to black 
holes with horizons. The equation for the horizon radius rn is given by solving the 
equation gu = 0. When we choose the positive root, gu can only vanish if the sum 
of the first two terms in (j3.7p is negative for values of r smaller than ru , in order to 
cancel the positive contribution of the square root as r tends to rn- This condition 
implies that 



Furthermore, we should mention that the above condition holds whenever c > and 
A; < due to an analogous reasoning, taking into account that we must take the 
solution with the negative root. 




(3.10) 



ck < 



rjj >-2{d- 3) K'^ck. 



(3.11) 
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We would like to emphasize the above point: once we set cfc < and choose the 
appropriate sign for the root, then the black hole solutions of our theory have a 
minimum horizon radius given by (j3.1ip . 



3.2 Vacuum Solutions 

If we consider the vacuum black hole solution, fi = Q = 0, it allows for horizons but 
no singularities whenever A; < 0. The horizon radius of these geometries is given by: 



2 1^1 



{d-2) 



2\D\ 2k2D 



(3.12) 



When A; > there is no horizon and the vacuum solutions correspond to regular 
anti-de Sitter space with spherical spatial slices. When A; < then the = Q = is 
not actually the vacuum solution, as there are negative energy black holes. 

3.3 Horizon Geometries 

Taking k = ^{d — 2) or k = Q gives rise to three qualitatively different asymptotic 
spacetimes. In particular, the asymptotic structure we obtain is as follows: 

k = +{d-2) ^ S'^ xm} (3.13) 
A: = ^ M"^ X (3.14) 

k = -{d-2) xM}, (3.15) 

where the M} denotes the timelike direction and denotes the spatial d-dimensional 
hyperboloid. We will be studying the thermodynamics of all three branches in the 
sections that follow. 



3.4 Horizon, Temperature and Extremality 

We will proceed to study the horizon structure and Hawking temperature Tfj of our 
black hole solutions. The horizon radius rn is given by the value of r that leads to 
the vanishing of the gu component of the metric. This condition implies that rn 
obeyj^: 

- At^^^^'"' - c ^^~]}^'^~^^ k\H'-' + ^rn' + / {th) - A^ = 0. (3.16) 

*It may seem that p.l6[) has solutions for fcc < that violate the bound (|3.1ip : however, these 
solutions correspond to the metric (|3.7p with the opposite root. 
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The above equation will have multiple solutions for rn corresponding to multiple 
horizons. The true event horizon is given by the largest positive root of (j3.16p . 

The Hawking temperature is most readily obtained by the well known result [H] 
which reads AttTh = (e~2^)'|r=-rj^- Applying this formula to our metric gives the 
following equation for the temperature: 



AttTh 



((i-l)riy2 2kD 
2 d- 1 



d-3 



d 



-kr 



H 



d-2 



4A 



+ -tTh'' + {d-4)fi + rH 
a 



d-3 



Ijr) 

rd-4 



0. (3.17) 



We can always choose to parameterize our thermodynamics with two of the four 
variables rn, Th, Q and fi, since we have a constraint equation relating the temper- 
ature to the derivative of gtt, as well as the above constraint that rjj must obey. The 
most convenient choice for our discussion will be th and Q. The reason for this choice 
is evident from the complicated form of I{r), which can only be written explicitly as 
a function of Q and ru- In this language, we can express the Th and fi as functions 
of ru and Q as follows: 



Th 



d-2 



4tt 



(d~l)rg2 2kD 
7? d-2 



(3.18) 



I d-1 
~^d-2 



kr 



d-2 



H 



V- + + /(..). (3.19) 



We can obtain the horizon radius for the extremal black hole by solving Th = 0. We 
will mostly be interested in the explicit form for the Q = Q extremal black hole with 
gtt component follows from ()3.7p : 



9tt{r) 



-2A 



k (d-l)r2 



d-2 



+ 



K2^A 



4L>2 



and horizon radius satisfying: 

gtt{re) = 

The parameter //g is obtained by evaluating (j3.19p at re with Q = 0: 

A , 



1 d - 1 



n'^d-2 



(d-l)(d-3) , 

1^2 + d"^^ 



(3.20) 
(3.21) 

(3.22) 
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3.5 Asymptotic Solution 

For far away observers, we can take the r — > oo limit of our metric to obtain an 
expansion in powers of r: 

e =lrI + ^-Z^,+Z±f^,+0[-j^A. (3.23) 



£^ d — 2 7-ci-2 j.2d-A ' I j,Ad-& 

where we have made the following relations between the asymptotic parameters and 
those in our action: 

did-l) dD , 

IJ, = -mK, (3.25) 
Q2 = 2{d-2)Kq^/g^. (3.26) 

We immediately note that the above asymptotic form is exactly that of the charged 
black hole with a negative cosmological constant. It is important, however, to recog- 
nize that the parameters of our theory are Q, /i and A and not g, m and and it is 
with these parameters that we will continue to work for the rest of the paper. 

Finally, we can write down the mass M and charge Q of the black hole as measured 
by a faraway observer: 



M = -S.xi^, (3.27) 



= El,xlM5l!. ,3.28) 

where we have that is the {d — l)-dimensional volume of a hypersurface with 
curvature given by {d — l)k. Note that c is not required to be small in the above 
expressions. In section we will compute the thermodynamic mass and charge of 
the black hole in the grand canonical ensemble and find that they agree favorably 
with the above expressions when |c| <^ 1. 
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PART II - Gauss-Bonnet Thermodynamics 



4 Global Stability 

Having completed the discussion of our black hole solutions, we can proceed to ex- 
amine their thermodynamic phase structure. Hawking and Page demonstrated in 
their seminal paper [15] that asymptotically anti-de Sitter Schwarzschild black holes 
are thermally favored when their temperature is sufficiently high, whereas the pure 
anti-de Sitter background was thermally favored for low temperatures. As the tem- 
perature is decreased there is a first order phase transition leading the black hole 
spacetime to the pure anti-de Sitter spacetime. Subsequently, the effect has been 
studied for various other asymptotically anti-de Sitter black holes, including those 
with Gauss-Bonnet terms [TOl [III Eg [Z3 ESI 123 EOl ED E3 IMl ES] • Our aim in this 
work is to contribute to this discussion the case for charged black holes, in theories 
with both Gauss-Bonnet terms as well as terms. We begin by reviewing and ex- 
panding the £ = case which was considered in [27[|3Tj. Our results are formulated 
in the language of the rn, Q parametrization, which is best suited for the subsequent 
e / analysis. 

4.1 Free Energy 

The free energy for our black hole solution is obtained from the thermal partition 
function, Z, of our theory. The Euclidean metric is obtained by continuing to imagi- 
nary time and identifying it periodically. The period is given by the inverse Hawking 
temperature. 

Using the saddle point approximation we can evaluate the partition function by 
evaluating the classical action of the black hole metric. In particular, the free energy 
is given as follows: 

F = -knTH log Z « -kBTnlE". (4.1) 

We should pause for a moment and note that we are working in the grand canonical 
ensemble, where we keep that electric potential fixed rather than the electric charge 
and thus the free energy we are computing is in fact the Gibbs free energy. We hope 
to study the thermodynamics of the canonical ensemble, where the electric charge of 
the black hole is kept fixed, in a future work. 

The Euclidean action evaluated for the black hole solution is in fact divergent. In 
order to regulate this divergence, we cut off our space at finite yet large r = r^ax 
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and subtract off the Euclidean action of the reference background before taking the 
cutoff back to infinity. We note that one can obtain the free energy of the black hole 
through the method of counterterms (see for example [41^ B2t [l3l [HI B5] ) , which is 
particularly convenient if there is no natural reference background to subtract. 

The reference background for the k > is given by the pure thermal anti-de Sit- 
ter space obtained from (j3.7p with A; > and /i = Q = 0, along with a uniform 
external potential taking the value of (minus) the internal electric potential at the 
horizon of the black hole [36]. Such a uniform electric potential leads to a vanishing 
Maxwell tensor, and consequently the Euclidean action remains that of the original 
thermal anti-de Sitter space. 

Our reference background for the case of k = is given by the pure thermal anti-de 
Sitter space obtained from (j3.7p with A; = and fi = Q = 0. Once again, there is a 
uniform external potential as in the k > casej^ 

Finally, for the case of black holes with curvature k < the reference background 
is rather different. One possible regularization procedure is that of (SlEHlllS], where 
we subtract the p.7|) solution with A; < and /i = Q = 0. This solution is still 
asymptotically anti-de Sitter space; however, it also contains a horizon and thus a 
fixed temperature. In order to avoid the problem of our reference background having 
a fixed temperature, we will choose our reference background to be the neutral ex- 
tremal black hole (j3.20p which can have an arbitrary periodicity in Euclidean time. 
Note, however, that black holes with non-zero external potential cannot decay to 
this reference background, since the reference background acquires a non-vanishing 
Wilson loop around the Euclidean time for an arbitrary external potential. Thus we 
cannot discuss Hawking-Page transitions in the k < case, and we will only study 

^However, we should make an important remark about the k = regularization procedure. For 
the case of vanishing Gauss-Bonnet coefficient it was shown in [3^ that there exists a vacuum solution 
with energy lower than pure anti-de Sitter space, known as the anti-de Sitter soliton, that is also 
asymptotically anti-de Sitter space. This solution can be regarded as the reference background for 
the Ifawking-Page transition as long as we compactify one of the directions of the Ricci flat black hole 
in order to embed it asymptotically in the anti-de Sitter soliton background. The thermodynamics 
with the anti-de Sitter soliton vacuum were studied with no Gauss-Bonnet coefficient in [37] and 
in for finite Gauss-Bonnet coefficient but electrically neutral black holes. In this paper we will 
not consider the anti-de Sitter soliton vacuum as a reference, we will only explore the Hawking-Page 
structure for the unmodified asymptotically anti-de Sitter charged Gauss-Bonnet black holes with 
F"* corrections. It would be interesting to explore the thermodynamics of the anti-de Sitter soliton 
reference background for the theories we are considering. 



17 



the local phase structure with our free energy regulated by (j3.20p . We refer to [56] 
for more details (see also [17]). 



There is one more subtlety as to how we regulate the Euclidean action. We must 
ensure that the thermal circle of both the black hole solution and the reference back- 
ground are the same at the cutoff radius r-max in order to have the same geometry 
at r = Tmax- This is possible since, even though the Euclidean black hole can only 
have a particular temperature to avoid conical singularities, Euclidean anti-de Sitter 
space or the Euclidean extremal black hole can have any temperature and no conical 
singularities. 

We are now ready to compute the free energy. It is very helpful to take the trace 
of (j2.4p which yields the expression: 



^^^R-{d+l)A + c{d-Z) 



(rf-3) 
4^2 



Ff^.F'^" + ci (d - 7) {F^.F'^n" + C2{d- 7) F^.F^^FxpFP^' = 0. (4.2 



Using the above result allows us to decompose our Euclidean action into the following 
two pieces: 

Itot = h + l2, (4.3) 



where we have defined: 
2 



h = 
h = 



(d-3) 
4 



R-2k 



-PhFi, 



(4.4) 



d-3 

i-i 



d'^+^x^/^ 



ci {F^.^F^^'^r + csF^.F^^Fa^F"^ « -f3HF2 (4.5) 



and Ph = . Notice that we have eliminated the Gauss-Bonnet term from our 
action. 

We regularize our Gibbs free energy by subtracting away that of the pure thermal 
anti-de Sitter space or any other appropriate reference background with the thermal 
circles identified at the cutoff boundary. This leaves us with the expression: 



F\ 



— lim 

(a — 3) Tmax^OO 



max / "I 

drr'^-^ I —R - 2A 



A(T-max) + Ao(rmax) 



1 



drr"^-^ ^Rn-2A 



(4.6) 
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The factors of /3h have canceled since they also appear as part of the Euclidean time 
integral. Rq is the Ricci scalar of the reference background solution. The vacuum 
horizon radius, r^ac, vanishes for A; > 0, whereas it equals r^. in ()3.2ip for A; < 0. We 
can now use the following relation for the curvature scalar: 



R 



1 



+ 



{d-l)k 



(4.7) 



to obtain 



-^1 Fvac 

2 



d- 3 



1 



where F,,, 



for /c > and 



d-2 



H 



a 



d- 3 



1 d-l 



d-2 



o y^^i vac 

2 a 



(4.9) 



for A; < 0. The Hawking temperature Th as a function of rn and Q is given in (j3.18p . 
We finally note that F2 = for the case of vanishing e in accordance with [31j, so 
that in fact (j4.8p is the total Gibbs free energy of our thermodynamic system. 

In the grand canonical ensemble we are considering systems with a fixed external 
electric potential, 4>. Imposing that (p be fixed to cancel the electric potential at the 
horizon of the black hole, (f)BH, leads to the following expression: 



{rH,Q) 



gQ 



2{d- 



„d-2 



{rH,Q)- 



(4.10) 



This expression will be of use when computing the thermodynamic energy, entropy 
and charge of the black hole. Furthermore, we will study the global and local stability 
of the system along slices of constant (p. 



4.2 Thermodynamic Quantities 

The Gibbs free energy contains all the information required to compute these quan- 
tities which are given by the usual thermodynamic relations. In particular, the ther- 
modynamic energy (which we will also call mass) is given by: 



E 



dh 
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Sfc X 



{d - 1) kr^jf^ ^ g^Q^rjf^ _ r|A ^ (d - 1) (d - 3) cPr'^f^ 



((i-2)K2 ' 2{d-2) d {d-2) 
^kxKrH,Q), (4.11) 



where we recall that Ii = —[5hFi was given in (j4.8p and our electric potential (j) is 
given in (j4.10p . The thermodynamic charge and entropy are given by: 

S = /3h(^)^-/, =Etx (i^l|L. + 8„t(<i-l)rJ-M. (4.13) 

Our thermodynamic electric charge has a factor of 2g that may differ from other 
expressions in the literature. This is due to our normalization of the kinetic term 
of the field strength in (j2.8p . We also notice that the thermodynamic entropy can 
become negative when cfc < (see for example |48] ) . In this case we find that the 
minimum radius of a black hole with positive entropy is given by: 

r1j{S = Q) = -2ck{d-l)K'^. (4.14) 

It is worth noting that since our free energy is parameterized by rn and Q, the 
above expressions are formal and should be computed as derivatives in rn and Q. 
We give the explicit expressions for the thermodynamic derivatives in appendix [Cl 
It is also a point of interest to notice that the thermodynamic energy and charge 
differ from the ones we obtain from the asymptotic metric in (j3.28p when c is not 
small. This should not come as a surprise since the faraway observables are not 
the conserved ADM charges of the metric. On the other, hand the thermodynamic 
quantities computed above do in fact agree with the conserved ADM charges (see for 
example [l9]) as expected. 

4.3 Hawking-Page Transitions 

Since we are parameterizing our thermodynamic variables with rn and Q, we must 
express ^ and Th in terms of the horizon radius and charge from (j3.17p by setting 
e = 0: 

^ K2((i-2) " {d-2) " d ^ 2((i- 2)4-2'^ ' 
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As in the usual Hawking-Page transition we are essentially going to search for 
regions in our parameter space where the free energy is negative, and identify these 
regions with black holes that are thermally favored over the reference background. 
We further note, as discussed in [10\ [TT\ I46j . that there may be regions where black 
holes are favored globally over the reference background, but are locally unstable 
in the sense that the specific heat of the thermodynamic system is negative for the 
globally stable configuration. This implies that even though we have located regions 
where the black holes are globally favored, they may still decay into some (unknown) 
configuration that presumably escapes our ansatz |10[ lllj . We will not attempt to 
find what these solutions are but we will still attempt to identify such regions in the 
analysis that follows. 

We should also mention once again that the e = analysis has been performed in 
\27\ [3T] , although we present it in slightly more generality since we are concentrating 
on arbitrary dimensions and arbitrary horizon curvature k. Furthermore, we present 
the results in the rjj, Q parametrization for d = 4 in order to readily compare the 
vanishing e results with the non-vanishing e ones. We have done the d > 4 analysis 
in the rn, Th parametrization and the results are presented in appendix [Dl 

In order to study the global thermodynamic stability of our system we can obtain 
an explicit form for the critical temperature for any c in any number of dimensions 
by solving (j4.8p at = for Th- When working in the th, Th parametrization for 
pure Gauss-Bonnet corrections, we show in appendix [D] that this results in: 



The significance of the critical temperature is as follows. For temperatures above the 
critical temperature with electric potential fixed at infinity, we find that the black 
hole solution is thermally favored globally with respect to the reference background 
solution with constant electric potential. Conversely, for temperatures below the 
critical temperature we will have that the reference background solution is globally 
favored. This is precisely the Hawking-Phase transition. 

We can also express the transition between globally stable black holes versus the 
reference background solution in terms of a critical charge which is relevant to the 



Tc = - 




(d-l)(d~3)fc^c 
d-2 



(4.17) 



47r 



kc {d — 1) rn 



21 



rH, Q parametrization we are working in. In particular, we find that the critical 
charge Qc is given by: 



Qc 



2{d-l) 



rkrH^ + cfe2 (^{d - 7) + 2c (d - 1) {d - 3) ^) 



l-2c(d-l) 



+ 



{d-2) Kth^ (l + Qc{d-l) 



d{d-l) (l-2c(d-l)^) 



2(d-3) 



(4.18) 



Whenever we are considering a configuration with Q > Qc we find that the black 
hole is globally favored thermodynamically, whereas if Q < Qc we have that reference 
background solution is globally favored. 

There arc two more critical values of the charge that are of interest. These 
correspond to those black hole configurations with zero mass and those black hole 
configurations with vanishing temperature. The zero mass black holes have a charge 
that satisfies the following expression: 



Q 



2 

M=0 — 72 



2 



id-1] 



krn^ - c{d-l)id-3) fc^ + 



{d-2 )A^ 



-rH 



2(d-3) 



(4.19) 



This bound is not physical in the sense that anti-de Sitter space allows for black holes 
with negative thermodynamic mass but positive horizon radius and entropy. However, 
we find it instructive to include this bound in our analysis of the thermodynamics, 
in order to explore regions of the parameter space containing negative and positive 
mass black holes with greater clarity. 

The second bound, which is in fact physical, is the cxtremality bound which 
requires that our black holes have a non-negative temperature. For a given rn the 
charge of an extremal black hole is given by: 



Qth=o 



9^ 



id - 1) ^ Jd-l)(d-3)(d-4)^, _ 



d-2 



rH 



2((i-3) 



(4.20) 



So the physical bound on the black hole charges following from the above expression 
is that Q < Qth=o- Before we study the full global phase structure expressed in 
the (f), c plane, it will be convenient to consider the local structure first and build 
the global picture from there. In the next section wc discuss the conditions for local 
thermodynamic and electric stability and study the various configurations arising. 
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5 Local Stability 



Having obtained the conditions for the global thermodynamic stability of a solution, 
we proceed to explore the local thermodynamic structure by computing the specific 
heat and electrical permittivity of our system. The specific heat informs us about the 
thermal stability of the black hole under temperature fluctuations, and the electrical 
permittivity informs us about the thermal stability of the black hole under electrical 
fluctuations. 

5.1 Specific Heat and Electrical Permittivity 

Even though a black hole configuration may be found to be in global thermodynamic 
preference to the reference background, it may still be locally unstable to some other 
globally favored configuration. It is not clear to what stable configurations such 
locally unstable black holes decay and our guess is that such configurations do not 
satisfy our simple ansatz. Additionally, we will also find configurations of multiple 
locally stable black hole solutions corresponding to a single temperature, for a fixed 
potential. Such configurations open the possibility for first order phase transitions 
from the metastable black hole to the globally favored one. 

In order to analyze local stability, we will consider two notions of local stability 
as discussed in ^Oj. Firstly, we will consider the specific heat of our system in the 
grand canonical ensemble. This is given by the following expression: 



For regions of the parameter space where the specific heat is positive, we have that the 
black holes are locally stable to thermal fluctuations. If the specific heat is positive, 
then an increase in temperature for fixed potential will result in an increase in the 
entropy resulting in a thermally stable situation. On the other hand, if temperature 
fluctuations lead to a decrease in the entropy of the black hole the system is locally 
unstable. 

Suppose for a moment that we have a temperature with n black hole solutions. 
Such a temperature will also have (n — 1) turning points as a function of the radius 
for fixed (f). For large radii we have that the temperature and entropy are both 
increasing functions of rn and thus the largest black hole has dS/drn > and 
{dTfj /drn)^^ > for a fixed (j). This condition leads to a positive specific heat 
and thus local stability under temperature fluctuations. Furthermore, the dTn/drn 




(5.1) 
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derivative changes sign consecutively with decreasing horizon radius, which in turn 
leads to a sequence of locally stable and locally unstable black holes for decreasing 
radius, so long as dS/drn > holds for fixed (j). The only case where we might have 
dS/drn < is for ck < 0; however, it is easy to show that such black holes will have 
a radius that violates the bound obtained in (|3.1ip and are thus rendered unphysical. 
We notice this feature throughout our analysis in the situations where we have many 
black holes. Finally, we note that if our temperature has only one black hole solution 
with positive entropy then this black hole will be locally stable. 

There is a second notion of local stability that we will study throughout the paper. 
This was introduced by Chamblin et al. in [lOj for the case of charged black hole 
thermodynamics and it is known as the isothermal electrical permittivity ex of the 
black hole. It is given as follows. 



dQ\ 



(5.2) 



If Et is non-negative our system is stable under electrical fluctuations. More 
specifically, if the electric potential at the surface of the black hole increases as a 
result of introducing an infinitesimal amount of charge to, the black hole the system 
is locally electrically stable. In contrast a black hole whose electric potential decreases 
when absorbing a small amoun. of charge is locally nnstablfl 

Having discussed our notions of thermodynamic stability, we can obtain their 
explicit expressions as functions of rn and Q. The specific heat for constant electric 
potential is given by: 



47r(d- l)r 
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(5.3) 



Let us clarify this with an analogy from classical thermodynamics. When studying the Gibbs 
free energy in a hydrodynamic system, our electric potential corresponds to the volume V of the 
system and the charge of the black hole to the pressure P. Then a positive isothermal electrical 
permittivity will correspond to minus the isothermal compressibility. In particular, a negative et 
would be analogous to a thermodynamic system that decreases in size as a result of increasing the 
pressure, which is clearly locally unstable. 
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The isothermal electrical permittivity is given by: 



2(d-2) (rj^ + k 




£t g 



1 + 
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We have defined the fohowing quantities in the above expressions: 



A ^ (d - 1) - (^2r|, + 
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(5.5) 



(5.6) 



At this point we can proceed to study the local stability of our black hole solutions 
for the various regions of our parameter space. For the sake of completeness we 
have also included the corresponding figures of the following analysis in the rn, Th 
parametrization for (i = 4 in appendix |Dj 

At this point, we should emphasize that the effects we uncover are qualitatively 
different from those obtain in the pure Einstein-Maxwell theory |101 [TT] . In other 
words, effects due to one-loop corrections become as significant as tree level effects, 
which suggests higher order effects are also important. It was shown in \2Q\ I27j . 
however, that for d = 4 spatial dimensions the Gauss-Bonnet corrections can be 
matched with gauge coupling corrections in a phenomenological dual matrix model. 
This gives some evidence that the thermodynamic phase structure we uncover should 
be visible in a similarly constructed phenomenological dual matrix model. 

5.2 /c> 

d = A Spatial Dimensions 

We begin by studying the c = in five spacetime dimensions, which was discussed 
in |101 I28j . When cj)^ < 3k /Sk^, we find temperatures corresponding to two black 
hole solutions. The smaller of these always has < and is thus rendered locally 
unstable under temperature fluctuations whereas the larger one is thermally and 
electrical locally stable. Only a subset of the thermally unstable are also electrically 
unstable; however, all thermally stable black holes are also electrically stable. When 
(j)'^ > 3k /Sk^ only one black hole solution exists which is always locally stable both 
electrically and thermally as well as being always globally stable. 

The case c > with d = 4 is depicted in flgures [2^a) and (b) for c < 1/6k^|A| 
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and c > 1/6k^|A| respectively [27]. By looking at the slices of constant (/) in figure 
[2]^a) we find slices which intersect the constant temperature slices once (for higher 
values of 0) or three times (for lower values of (p) corresponding to one or three black 
hole solutions. In the region with three black holes, only the largest and smallest ones 
are locally stable, whereas the intermediate sized black hole always resides within the 
region of local instability enclosed by the C^-curve. These black holes may also live 
in the region of electrical instability which resides within the region of local instabil- 
ity for all the cases we consider. The smallest black hole lives within the region of 
global instability enclosed by the global stability (F) curve and will thus always have 
a decay channel to thermal anti-de Sitter space. The largest black hole may or may 
not be globally favored. A configuration where a black hole is globally unstable may 
have two decay channels: either the reference background or the black hole with lower 
free energy. Such metastable configurations reside within the region enclosed by the 
metastable (MS) curve. For lower values of the potential, the constant (j) curves cross 
the -F-curve, thus it can experience a Hawking-Page transition. For larger values of 
(j), the constant slices cross the extremality curve, Th = 0, without ever reaching 
the i^-curve. This means, in particular, that even at zero temperature the ground 
state of the boundary theory has non-vanishing thermodynamic entropj0. 




As we increase the value of c, we no longer find temperatures with three black 

^Note, however, as discussed in [21] that this does not imply that there is no longer a confinement- 
deconfinement phase transition. The temporal Wilson lines have vanishing expectation value since 
they are unable to wrap around a thermal circle with infinite periodicity. Similarly the spatial Wilson 
loops are not obstructed by the horizon which is infinitely far down the throat for an extremal black 
hole and thus cannot exhibit an area law. 
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hole solutions over slices of constant (p. Instead we obtain the situation in figure [2jb). 
The region of local instability has disappeared altogether and we find single black 
hole solutions that can undergo Hawking-Page transitions. Furthermore, for large 
enough <j) some of the black holes in figure [2)[^b) are globally favored along the whole 
constant (j) slice and experience no Hawking-Page transition. 

Finally, we consider the case of c < for d = A. When c < -l/6|A|ft;^ the be- 
havior is identical to the d > 4 case and it will be described below. Here we describe 
the case c > — 1/6|A|k^, which is depicted in figure [3l For slices of constant (f) we 
find constant temperature slices that intersect them at either one point for larger 
values of (j) or at two points for smaller values of (j). For regions with two black holes, 
the smaller black hole is always within the region of local instability. There are two 
cases to be noted. The first case corresponds to slices of constant 0, where the larger 
stable black hole can experience a Hawking-Page transition by moving along the (p 
slice. The second case corresponds to the larger black hole unable to experience a 
Hawking-Page transition because it encounters the region of negative entropy, to the 
left of the 5'-curve or C^-curve, before reaching the F-curve. The latter case occurs 
for relatively large values of (j). Finally, we find that for those slices of constant tem- 
perature intersecting the constant (/> slices at a single point, and thus corresponding 
to a single black hole solution, there is no Hawking-Page transition accessible. More 
specifically the constant cf) curve ends at the extremality bound without intersecting 
the global stability bound and thus remains globally stable as we continuously lower 
the temperature. 
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Fig. 3: Thermal structure for k = 2, d = A and c 
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d > 4: Spatial Dimensions 



For the higher dimensional case we find that the c = case remains unchanged 
^LU[ [28] . Once again there is a critical value of the electric potential (j)^ = 2(d-2)^K^ • 
If (/)^ > we have a black hole solution that is always globally and locally stable 
and thus undergoes no Hawking-Page transition. If (fp' < cf^ we have two black hole 
solutions of which only the smaller one is locally unstable and the other may undergo 
Hawking-Page transitions as we vary the temperature. 



When c > in higher dimensions, we find temperatures with two black hole solu- 
tions of which the smaller one is always locally unstable. For small enough values of 
the potential we have that the larger black hole may undergo a Hawking-Page tran- 
sition. For large values of the potential the constant (p slices end at the extremality 
bound, as is evident from figure H] and thus there is no Hawking-Page transition. 
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Fig. 4: Thermal structure for k = d — 2, d > 4 and c > 0. 



Finally, when c < in higher dimensions we find various deviations from the case 
with five spacetime dimensions. Let us begin with small negative values of c as ex- 
hibited in figure [5](a) . We find temperatures with three black hole solutions of which 
the smallest and largest are always locally stable, whereas the intermediate sized 
black hole is always locally unstable. For such three black hole configurations the 
smallest is always globally favored, but we find values of (p where all three can be 
globally favored for a particular temperature. At such temperatures we can posit 
the existence of a metastable phase transition to the black hole with the lowest free 
energy provided that it has positive entropy. In the case that the globally stable 
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black hole has negative entropy such transitions cannot occuijj. This phenomenon 
happens for arbitrarily small negative values of the Gauss-Bonnet coefficient c and 
should in principle be visible in the thermodynamic structure of the dual boundary 
theory. Once again the region of metastability is bounded by the MS'-curve and the 
C^-curve and we find that either the larger or smaller black holes can be metastable. 
Furthermore, the large metastable black holes that live within the metastable region 
but above the MS, 5-curve would decay to small black holes with negative entropy. 
Thus such large black holes are not truly metastable modes. 

It is perhaps amusing that we also find slices of constant (j) that cross the global 
stability curve at two points. Naively this would seem like the possibility for two 
Hawking-Page transitions. It turns out that only one of the Hawking-Page transitions 
occurs outside the region of local instability; only the largest black hole configuration 
can experience a Hawking-Page transition as the temperature is varied. We also find 
constant temperature slices with two black hole solutions. Only the larger one is 
locally stable and can always experience a Hawking-Page transition upon varying the 
temperature. For large values of 4> we find single black hole solutions that are always 
globally and locally stable and never cross the i^-curve as can be seen from[S|^a). 




(a) (b) 

Fig. 5: Thermal structure for k = d — 2, d > 4. The two regions of c displayed are: 

(a) c < c < and (b) - ^a%fd%tlf-5)%A\ < c < c. 

As we push c to slightly more negative values, we find that for any value of (j) with 
the three black holes these are always globally stable and thus undergo no Hawking- 
Page transition. Furthermore, there are regions where the smallest of the three black 
holes is in a region of negative entropy so that the metastable phase transition towards 

*See [31] for a different interpretation about tfie negative entropy black lioles 
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these cannot occur. These configurations are displayed in figure [5|^b). 

As c is cranked down to even lower values we find that the two smaller black 
holes occur in regions of negative entropy and the metastable transitions are lost 
altogether as is observed in[6ja). We also find slices of constant (j) intersected twice 
by the same temperature, corresponding to two black hole solutions. A small subset 
of these slices of constant (j) encounter the global stability curve at three points. 
Once again, only the largest black hole configuration is physical and can experience 
a Hawking-Page transition. The other two Hawking-Page transitions would happen 
in regions where the black hole is either locally unstable or has a negative entropy 
and hence are rendered unphysical. These features are shown in figure [6l^a) where 
the two unphysical Hawking-Page transitions occur in a region of local instability 
and figure [Hb) where one of the unphysical Hawking-Page transitions occurs in a 
region of negative entropy yet local stability and the other unphysical Hawking-Page 
transition occurs in a region of local instability. 




(a) (b) 



Fig. 6: Thermal structure for k = d — 2, d > 4. The two regions of c displayed are: 

(d-3)(d-4)(2d-3) d{d-3) , 

^"^^ 4((i-l)(d-2)(2d-5)K4|A| ^ t- ^ 4(d-l)(d-2)K4|A| ^'^'^ 

(U\ d(d-3) ^ r ^ 

^^1 4(d-l)(d-2)K4|A| ^ ^ ^ 4{d-2)K4|A| • 

We encounter values of c where there are either one or two black hole solutions as 
seen in figuresEl^b) and[71^a) and[7|^b). In figure[7[a) we find that the constant slices 
with two black holes can cross the global stability curve three times for larger values 
of (/> or once for smaller values of (/>. Once again only the Hawking-Page transition 
of the largest black hole is physical. In figure El^b) there is no unphysical Hawking- 
Page transition for constant ^ slices with two black holes. The constant (\) slices with 
single black hole solutions either cross the global stability curve twice for smaller 
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V-2(d-3)ck K 




V-2(d-3)ck K 4.5 



(b) 



Fig. 7: Thermal structure for k = d — 2, d>A The two regions of c displayed are: 

d{d~l){2d~5) d „^ J (Y,) „ ^ d{d-l){2d~5) 

^"^^ 4(d-2){d~3)(2(i-3)K4|A| ^ ^ ^ 4(d-2)K4|A| ^^^^^ ^ 4.{d-2){d-3){2d~3)K^\A\ ' 



values of or never for larger values of (p. In the case that they cross it twice, only 
the Hawking-Page transition between the reference background and the largest black 
hole configuration is physical. 

5.3 k<{) 

We now turn our attention to those black holes with hyperbolic horizons. The sit- 
uation with c = was studied in [38], where it was found that the black holes are 
both globally and locally favored and experience no Hawking-Page phase transition 
in the grand canonical ensemble, as we discussed in section WA\ This result holds for 
aU d > 4. 

The case of c > is depicted in figures [8] and [9l The reference background, which we 
have chosen to be the neutral extremal black hole, resides at the point RB. The first 
situation we encounter is shown in figure [8l^a), where we have temperatures with sin- 
gle black hole solutions. Such black holes can have negative or zero mass only if they 
reside within the region bounded by the M-curve. The next situation which is closely 
related is given figure [5]^b) . In this case we still have temperatures corresponding to 
single black hole solutions; however, we also encounter the S'-curve discriminating 
those solutions with negative and positive entropy. This entropy bound, however, 
does not enclose all the massless and negative mass modes. 

As we increase the Gauss-Bonnet coefficient c we find that the negative entropy 
bound has engulfed all the massless and negative mass modes, as seen in figure [9)[a). 
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\/-2(d-3)ckir V-2(rf-l)ck«: 



(a) 



V-2(d-3)ckir V-2(d-l)ckir 

(b) 



Fig. 8: Thermal structure for k = — {d — 2), d > 4. The two regions of c displayed 
die. i^a; c ^ 4((i-i)(d-2)Ki|A| '^'^'^ '^"^ 4.{d-l){d-2)K*\A\ ^ ^ 4{d-l)(d-2)ft4|A| • 



We also arrive at a similar situation in figure [9)[|b), where we have temperatures with 
two black hole solutions as well as temperatures with single black hole solutions. 
As is usual the smaller of the two black holes is always locally unstable as well as 
having negative entropy, whereas the larger one is always locally stable. There are 
no massless modes present with or without negative entropy in this situatior§|. 
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Fig. 9: Thermal structure for k 

(^) 4{d-l){d^2)K^\A\ < C < 



— {d — 2), d > 4. The two regions of c displayed 
djd-l) J /i^N „ . djd-i) 



®One may be interested in possible decays between black holes and the charged extremal black 
hole with the same external potential. We find that such decays are only possible when the extremal 
black hole possesses negative entropy. 
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For the case c < we find that there is again a critical value of (j) that separates the 
phase structure into two regions. The region with small values of (j) has temperatures 
corresponding to black hole solutions that can have zero, negative or positive mass 
and the region with large values of (j) where the black holes have positive mass. This 
situation is qualitatively equivalent to the one in figure [8]^a). 

6 Global Thermal Phase Structure for c ^ 8z e = 

6.1 k>0 

d = 4 Spatial Dimensions 

We begin by summarizing the thermal phase structure in the (j) vs. c plane for 
d = 4 spatial dimensions. The situation is depicted in figure [TUlfa) and a derivation 
of the curves bounding the various regions is provided in appendix |El There are 
three main regionj^ denoted by the Roman numerals I, II and /// containing one, 
two and three black holes respectively. These regions are identified by mapping out 
the regions of local instability given by the region within the green curve. Once 
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Fig. 10: Thermal phase structure for A; > and d = 4. 
we have identified the three main regions we can identify subregions given by the 



'For convenience we have also included the definitions of the various regions in Table 1. 

33 



particular properties of the configurations. These are labeled by the subscripts and 
superscripts. The subscript E denotes regions where there are extremal black hole 
configurations. The subscript denotes regions where the only extremal solutions 
are ones with vanishing horizon radius. The subscript S denote regions where there 
exist configurations with negative entropy. In S regions, possible extremal black holes 
always have negative entropy. Finally, the superscript HP denotes configurations 
that experience a Hawking-Page transition. The superscript HP* denotes the region 
where the Hawking-Page transition for the large black hole happens in a region of 
negative entropy. Note that we have metastable black holes in region /// where there 
are three black hole solutions. Of particular interest is region /// where we have three 
black hole solutions of which the intermediate one is locally unstable. The metastable 
configuration can decay to either the more stable black hole hole configuration or the 
reference background. 

d > 4 Spatial Dimensions 

We now proceed to the d > 4 case depicted in figure [TTJ Once again, the regions of 
the diagram are obtained explicitly in appendix[El The region within the shaded box 
in figure [TT] is expanded. Once again the phase diagram is divided into three main 
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Fig. 11: Thermal phase structure for k > and d > 4. 

regions /, // and /// denoting the number of black hole solutions. These regions 
are delineated by the curves of local instability. As before the metastable black holes 
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are found in region Within each main region there are various subregions as 

before. We have encountered, however, two new types of subregions. The subscript 
S* denotes a subregion of /// where the smahest black holes always have negative 
entropy. The subscript HP** denotes subregions where we find the lines of constant (f) 
cross the global stability three times and thus there are two unphysical Hawking-Page 
transitions. Of particular interest is the region IIIe^Si where we encounter three 
globally stable black hole solutions at a constant temperature. The intermediate 
sized one is locally unstable. If we are in the metastable configuration we can have a 
phase transition to the globally stable configuration without decaying to the reference 
background. 

6.2 k<0 

We proceed to study the case with k < which is depicted in figure [T2l As before 
the explicit regions given in these diagrams are derived in appendix [E) There are 
only two main regions denoted by / and //, where we have one or two black hole 
solutions. As usual, these regions are split into various subregions. A new subregion 
that we have not yet encountered is denoted by the subscripts M and M*, where 
we find black holes with vanishing or negative thermodynamic mass. The asterisk in 
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12: Thermal phase structure for < 0. 
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M* refers to regions where the negative or vanishing mass black hole has negative 
entropy. Finally, we mention that black holes that are unstable locally occur only in 
regions of parameter space that are hidden by the negative entropy bound. 

6.3 k = 

Finally, we visit the k = case which has the simplest thermal structure. Particularly, 
for the case of a flat asymptotic geometry in the grand canonical ensemble, we find 
again the absence of a Hawking-Page phase transition. This is the same result as 
that for the c = case and we find that it holds for all values of c. 



Table 1: Definitions of Subregions 


/ 

// 
/// 


Single black hole for any given temperature. 
Two black holes for any given temperature. 
Three black holes for some temperatures. 


HP 

* 

** 


Physical Hawking-Page transition. 

Single unphysical Hawking-Page transition. 

Two unphysical Hawking-Page transitions. 


E 


Extremal black hole present. 

Extremal black hole has vanishing radius. 


S 
S* 


Negative entropy modes present. 

Smallest black holes have negative entropy. 


M 

M* 


Negative mass modes present. 

All negative mass modes have negative entropy. 



As a final remark we reiterate here that there is a vastly richer thermal phase 
structure for the case of charged black holes with a Gauss-Bonnet term than the 
usual Hawking-Page transition between neutral black holes and thermal anti-de Sitter 
space. There are configurations with locally unstable black holes which are globally 
stable as those discovered in [10l[Tl], as well as globally unstable yet locally stable 
configurations. The black holes we are considering have one more free parameter 
compared to the neutral black holes considered in the original work [15]. Thus the 
free energy is a function of two varying parameters constrained along slices of constant 
(j). A metastable black hole can thermally evolve to some r|^, Q* configuration that is 
locally and globally stable. This is particularly true in the grand canonical ensemble, 
where the electric potential is kept fixed but not the electric charge so that the black 
holes can emit and absorb charged radiation. 
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PART III - Thermodynamics 



7 Global Stability ior e ^ and c = 

We can now direct our attention to the case of e 7^ 0. Once we obtain the general 
expressions for the free energy, specific heats and thermodynamic quantities of the 
theory, we proceed to explore thermal phase structure of the theory with no Gauss- 
Bonnet term in the action. Finally, we study the deformations to the thermodynamic 
structure of the Einstein-Maxwell-Gauss-Bonnet theory caused by corrections. 

7.1 Computing the Free Energy 

In order to compute the free energy of the e 7^ theory we need to slightly extend 
the discussion in section 14.11 In particular, we have to evaluate the full Euclidean 
action given by the sum of (j4.4p and (j4.5p . Once again we will be working in the 
rH,Q parametrization due to the non-trivial Q dependence of the Euclidean action. 
The first term in the Euclidean action becomes: 



Ff - 

-'I vac 



d-3 



'•h 



(7.1) 



where now Th and ^ are considered as functions of rn and Q and furthermore depend 
on e as dictated by ()3.18p and (I3.19p . F^^^ was defined in ()4.9p . The second term in 
the Euclidean action becomes: 



d-3 



drr^~^f{r)\ (7.2) 



The above expressions are for general e and no approximations have been made. 
We will perform our analysis numerically using the numerical integration packages 
contained in Mathematica [50]. Thus the basic object we will be working with is 
given by the full Euclidean action: 

iLai = n + 12 = -m + Fi). (7.3) 

The reference vacua used to test the global stability of a given black hole solution are 
once again those introduced in section 14.11 since once we set Q to zero all e effects 
vanish. For the sake of completeness we include the x ^ 1 expansion of the total free 
energy below: 
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(7.4) 



Furthermore, we include the Hawking temperature expanded to first order in x- 
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(7.5) 



7.2 Thermodynamic Quantities 



Having obtained the Euchdean action we can in principle obtain all the thermo- 
dynamic quantities. We must compute the appropriate thermodynamic derivatives 
given in appendix[C]for the r//, Q parametrization. We note, however, that the fixed 
external electric potential is modified to the following: 



(7.6) 



which in the limit X ^ 1 becomes 
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(7.7) 



Even though we are barred from obtaining full analytic expressions, by direct appli- 
cation of Wald's formula [51j . we expect that the entropy as a function of rn and 
Q is not affected by a contribution to the matter Lagrangian. Indeed we can verify 
this to third order in x and numerically for arbitrary values of e. In fact one can 
explicitly check that the thermodynamic electric charge also remains unchanged for 
non-zero e. Thus we obtain: 
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We can then use the first law to obtain the thermodynamic energy 



total 
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(7.10) 
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Using the expansion in (j3.6p we expand /i^ in x ^ 1: 
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X + 0{x^) . (7.11) 



7.3 Hawking-Page Transitions 

In principle we would like to obtain the analjd;ic quantities for the critical parameters 
that dictate the Hawking-Page transitions. We can achieve this by expressing I^qIqi in 
terms of our two free parameters th and Q and solving /fo^^; = 0. This is only possible 
numerically for arbitrary values of e although one can obtain analytic expressions 
upon expanding in x, which are given below: 



where Ql, Q\j=o and Q^^^q are given in (I4.18p . (14.19P and (j4.20p respectively where 
we set c = 0. We have also defined the unitless parameter Xx = (S^^Q^^)/?'^'^ ^\ 
where x = {c,Th = 0, M = 0}. Black holes with Q > are thermodynamically 
favored, massless black holes have Q = Q%j^q and positive temperature black holes 
satisfy Q < Qt^^q- 

8 Local Stability 

To test the local thermodynamic stability for non-zero e we must obtain the specific 
heat, Cff), and the isothermal electrical permittivity, et, of our solutions. These 
quantities are defined by equations (|5.3p and (|5.4p . Analytic expressions for C^j, and 
Et can be obtained only upon expansion in e and we will treat them on the same 
numerical footing as everything else. For completeness, however, we include the 
expansions for the local equilibrium parameters below: 




(7.12) 



(7.13) 



(7.14) 
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X + 0{x'). (8.1) 



Having obtained all the basic global and local thermodynamic quantities, we can begin 
exploring the various regions of the parameter space. Our approach will consist of 
first exploring the thermal phase structure for the theory with only an correction 
and no Gauss-Bonnet term. Once we identify the relevant qualitative regions of the 
e — (/>^ plane with c = 0, we will proceed to explore how the term deforms the 
various qualitative regions we found in the previous section for the e = case. 



8.1 fc>0 

We begin by exploring the e > case, for d > 3, depicted in figure [T3l For large 
enough values of (j)^ , we find that there is only one black hole which is always glob- 
ally and locally favored and undergoes no Hawking-Page transition. As we decrease 
the value of (jf , we encounter constant (p^ slices with three black hole solutions. Of 
the three black holes the medium sized one is always locally unstable, whereas the 
largest one can experience a Hawking-Page transition. The smallest black hole is al- 




Fig. 13: Thermal structure for > 0, d > 3 and e > in both the ru, Q plane (a) 
and the th, Th plane (b). 

ways globally favored and it cannot reach the global stability curve via a continuous 
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change in the temperature without first crossing the local instability curve. Interest- 
ingly we find temperatures where two of the three black hole solutions are globally 
favored and separated by an intermediate black hole which is both locally and glob- 
ally unstable. For such situations we can conceive of a phase transition from the 
least to most globally stable configuration. We also find temperatures where all three 
black hole solutions are globally favored although as always the intermediate one is 
locally unstable. A similar phase transition should happen for such situations. The 
metastable black hole is qualitatively similar to the one we found in the c < 0, > 
case for d > 4 with e = 0. The local electrical stability curve is always hidden within 
the local thermal instability region, as was the case for the c 7^ 0, e = black hole. 
This is in fact a common feature for all values of e and k both positive and negative. 

For the e < case, with d>3, depicted in figure [TU^ a) and (b), we find constant (j)^ 
slices containing one or two black hole solutions. For large enough values of (p^ we find 
a single black hole solution that is globally favored and undergoes no Hawking-Page 
transition. Lowering the value of (p^ we find that the single black hole solution has 
access to Hawking-Page transitions. Finally, for small enough values of (f)^ we have 
two black hole solutions for which the larger one experiences Hawking-Page transi- 
tions and the smaller is always locally unstable. We find that all extremal black hole 
solutions disappear after a critical e < due to the bound on the minimum allowed 
horizon radius (j2.19p depicted by the e-bound-curve in figure [Hlfbl. This bound also 
reduces the region of allowed Hawking-Page transitions as e — > — 00. 
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8.2 k<0 

The case for e > 0, with d > 3, shown in figure [15] constitutes of a single black hole 
solution which may be extremal, that is always globally and locally stable and never 
undergoes a Hawking-Page transition. A similar situation holds for e < displayed 
in figures [TBlfa) and (b), except that there is a further constraint due to the bound 
(j2.19p that eliminates extremal solutions for large enough Finally, for the A; = 
case there is as usual no rich thermodynamic structure. 

Q 




Fig. 15: Thermal structure for k < 0, d > 3 and e > £2 • 




9 Global Thermal Phase Structure for £ 7^ &; c = 

Let us summarize the results for e 7^ and c = with the global thermal phase 
diagrams in the e — ((/»^)~^ plane. 
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9.1 k>0 

We begin by examining the case with k > depicted in figure [T71 Notice that we 
have expressed the thermal phase structure in the e — {(p^)~'^ plane, so that large 
potentials are in the lower part of the diagram. The phase diagram is separated into 
three main regions with one or three black hole solutions and denoted by /, // and /// 
respectively. The expressions for the curves bounding the various regions are provided 
in appendix|El Of particular interest is region /// where we encounter the possibility 
of three globally favored black holes leading to phase transitions between the least 
thermally favored to the most thermally favored. This situation is qualitatively the 
same as the one we studied in the c<0, e = 0, A;>0 case for d > 4 spatial 
dimensions. However, it is a new phenomenon in d = 4 spatial dimensions, since in 
the c>0, e = 0, A;>0 case the smallest black hole always has a higher free energy 
than the reference background. 

1 




Fig. 17: Thermal phase structure for A: > 0. 

9.2 k<0 

The case with A: < given in figure [TH] is much simpler. There is only one main region 
containing single black holes and denoted by I. We find subregions with extremal 
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10 Metastability 
10.1 Gauss-Bonnet 

We now discuss the various metastable configurations we have encountered through- 
out the analysis. We begin with the e = case depicted in figure [T9l The results are 
expressed as constant (p slices in the F — Th plane. The locally unstable black hole is 
always depicted by the dashed curve connecting the small and large black hole curves. 
The d = A case, shown in figure [TW a) . clearly shows that one of the black holes is 
always globally unstable and corresponds to region III§o' in figure fTOTa). Thus the 
metastable black hole has two decay channels; the reference background or the ther- 
mally favored black hole. Figures [TWb). (c) and (d) are the various situations we 
encounter in the d > 4 case. The black holes have negative entropy to the left of the 
black dot. In figure [TWb) we find configurations with two decay channels as well as a 
single decay channel corresponding to regions Illj^g* and IH^s 5* of figure fTOTb). 
In figure [19] (c) and (d) we have a single decay channel since the metastable black 
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hole has lower free energy than the reference background. Thus these metastable 
configurations are found in regions IIIe.s and IIIe,s,S* of figure [TOTb) . 




(c) (d) 
Fig. 19: Metastability for c / 0, e = with (a) d = 4 and (b), (c) and (d) d > 4. 

10.2 

The c = situation is depicted in figure [20l We find configurations where the 
metastable black hole can decay to either the reference background or the favored 
black hole corresponding to region 1 11^^* in figure \T7\ or configurations where the 
metastable black hole has a single decay channel to the globally favored black hole 
corresponding to region I He in figure [T71 
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Fig. 20: Metastability for e / 0, c = 0. 

11 Deformations of Gauss-Bonnet Thermodynamics 

Our remaining task is to explore how the interesting phenomena we found in the 
thermal phase structure of the Gauss-Bonnet black hole are deformed by turning on 
the £ coefficient. Our analysis is by no means exhaustive of all possible cases; however, 
various new phenomena are uncovered when the e coupling is of comparable order to 
the c coupling. 

11.1 d = A Spatial Dimensions with e > and c > 

We begin with the c, e > case with four spatial dimensions. Upon deforming the 
regions of c with three black hole solutions by turning on the e, we find that tem- 
peratures for which the two locally stable modes can in fact become globally stable. 
This situation differs from the e = case where the smaller locally stable black hole 
was always thermally unfavored with respect to the thermal anti-de Sitter vacuum. 

We also find constant (j) slices that intersect the global stability curve at three 
points. Of the three intersection points, one is always within the region of local 
instability. Interestingly we find that for the other two points the Hawking-Page 
transition occurs in a physical region of parameter space. This differs from our 
previous encounter of such multiple Hawking-Page transitions, which always occurred 
in unphysical regions. We depict the situation in figures [2]Ta) and (b). Notice that we 
if we begin in a globally stable black hole configuration and continuously decrease the 
temperature the system will only ever experience a single Hawking-Page transition, 
since the two regions where Hawking-Page transitions are available are separated by 
a region of local instability. The constant (j) slices in the F — Th plane are observed 
in figure ED it is clear that for some slices both the large and small black hole can 
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(a) (b) 

Fig. 21: Thermal structure for A; = 2, d = 4 and c = 5Q^^|7q> £ = ^^^\x\ 
relevant metastability curves (b). 



experience a Hawking-Page transition. 

11.2 d > 4 Spatial Dimensions with e > and c > 

For the c > 0, e > case with five or more spatial dimensions, depicted in figure 
[22} we find temperatures intersecting the constant (j) slices four times, corresponding 
to four black hole solutions at a given temperature. Though novel, the smallest 
black hole has negative specific heat in accordance to our general argumentation. 
Furthermore, constant (f) slices that have two regions with physical Hawking-Page 




Fig. 22: Thermal structure for k = 3, d = 5 and c = gQ^4|^| , e = jQ^rpq- 

transitions. These regions are separated by a region of local instability so one cannot 
have two consecutive phase transitions as we continuously lower the temperature from 
the globally favored large black hole configuration. We see that phase transitions at 
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constant temperatures between two globally favored black hole states occur for c > 
in d > 5 dimensions with positive e switched on. Such phase transitions were only 
visible for c < in the pure Gauss-Bonnet case. 

11.3 d = A Spatial Dimensions with e > and c < 

The c < 0, e > situation for d = A spatial dimensions is given in figure [23l The 
situation is qualitatively identical to the c < case with d>b uncovered in the pure 
Gauss-Bonnet theory and shown in figures [5]^ a) and (b). There are constant (j) slices 
with three black hole solutions at a given temperature of which the smaller and larger 
may be globally favored over the thermal anti-de Sitter vacuum. We simply note that 
this situation had not been observed for the pure Gauss-Bonnet case with c < and 
d = 4. No secondary Hawking-Page transitions are observed in this case. 



Fig. 23: Thermal structure for k = 2, d = A and c = — 5q^4|^| , e = ^gi\{^\ • 

12 Discussion and Summary of Results 

We have uncovered various novel features in the thermal phase structure of both 
the Gauss-Bonnet black hole with no F'^ corrections and the Gauss-Bonnet- black 
hole. Our results are contained in the global phase structure figures [TOl [TT] and [12] 
for pure Gauss-Bonnet corrections, and in figures [T7] and [18] for pure corrections. 
In section [TT] we discus the case with both couplings turned on. We have mapped 
out various regions of physical interest within the phase diagrams, such as regions 
containing extremal black holes, regions including massless and negative mass modes, 
and regions with negative entropy modes. We depict all possible phase transitions 
between locally stable black holes and/or thermal anti-de Sitter space in figure [Ml 
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Pure Gauss-Bonnet 
d=4, oO 



LBH 
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^TAdS^HP 

Pure Gauss-Bonnet 
d>4, c<0 
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LBH 

HP 
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Pure F*, d>=3, e>0 



SBH 



LBH 



H^<^TAdS-^HP 
Gauss-Bonnet-F* 
d>=4, £>0, c<0 



Fig. 24: All possible physical phase transitions in the Gauss-Bonnet- theory. SBH 
denotes the small locally stable black hole, BBH denotes the large locally stable 
black hole, TAdS denotes thermal anti-de Sitter space, and HP denotes an accessible 
Hawking-Page transition. 



For the pure Gauss-Bonnet case we have discovered new metastable modes for k > 
and c < in d > 5 spatial dimensions, as depicted in figures [TWb). (c) and (d). 
More specifically, we find temperatures at a constant electric potential with three 
black hole solutions, of which the intermediate sized one is locally unstable. The 
smaller black hole is an 0{a') effect and thus its size is also of order The larger 
black hole has a size of the order of the anti-de Sitter length. These modes have the 
interesting property that for certain values of the electric potential, they may both 
be thermally favored over the thermal anti-de Sitter vacuum at a fixed temperature. 
This differs from the situation in d = 4 spacetime dimensions |27j . shown in figure 
[TWa). where one of the metastable modes is always thermally unfavored with respect 
to the thermal anti-de Sitter vacuum. It may be somewhat alarming that the effect 
occurs for c < 0; however, from a purely effective field theory point of view, we must 
take these effects seriously, unless there is evidence that this is an unphysical region 
of parameter space (see [52l [53] for a relevant discussion for the neutral case). It 
would be very interesting to understand the dual field theory interpretation of these 
metastable modes given that the phase transition occurs between two states in the 
deconfined phase. 

Another curious feature we uncover is the existence of constant electric potential 
slices that intersect the free energy curve more than once. Such a scenario would im- 
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ply that as one lowers the temperature, the large hot black hole decays into thermal 
anti-de Sitter space, and the thermal anti-de Sitter space itself eventually collapses 
into a black hole at even lower temperatures. This, however, is not the case, since 
the secondary Hawking-Page transitions always occur in either a region of local in- 
stability or a region of negative entropy, and are thus unphysical. 

For the pure F'^ case, which has not been previously explored, we uncover additional 
metastable modes for d > 3, k > and e > 0. Once again, they consist of three black 
hole solutions, of which the intermediate sized one is locally unstable. Furthermore, 
we find that for certain values of the potential, both locally stable modes are ther- 
mally favored over the thermal anti-de Sitter at a fixed temperature, as depicted in 
figure [20l and thus they are both states in the deconfined phase of the dual theory. 
This feature persists for d = 3 spatial dimensions when there is no Gauss-Bonnet 
correction present. Once again, any secondary Hawking-Page transitions observed 
occur in unphysical regions of parameter space. 

Finally, we study the thermodynamics of our theory with both couplings turned on in 
section [TTl For d = 4 and both couplings positive, we uncover constant (j) slices allow- 
ing for three black holes at a constant temperature, of which the smallest and largest 
may be globally favored. Furthermore, the secondary Hawking-Page transition be- 
comes physical; thus, for a slice of constant (p there are two disjoint Hawking-Page 
transitions that can occur - starting from either the small or large black hole. 

For d > 5 there are temperatures with four black hole solutions for a constant 
4>, indicating the existence of a fourth locally unstable black. Furthermore, the sec- 
ondary Hawking-Page transition for the small black hole is still present. It is of 
great interest to find a clear interpretation of the additional Hawking-Page transition 
from the point of view of the dual theory. Naively, it would seem that there are 
two disconnected deconfined phases at different temperatures that can experience a 
Hawking-Page transition to a confined phase as we continuously lower the tempera- 
ture. However, such phenomena deserve a complete treatment in their own right. 

We reiterate that a qualitative change in the tree level thermal phase structure arising 
from 0{a') corrections may imply that we must take into account all higher order 
corrections to get the correct physical picture. However, it has been possible to match 
the Gauss-Bonnet effects in d = 4 with the thermal phase structure of a phenomeno- 
logical dual matrix model |26t [27] with no corrections. This observation leads 
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us to expect that a similar situation occurs for the thermal phase structure we have 
uncovered. It would be worth further studying this issue by obtaining the thermal 
phase structure of effective gravitational actions with even higher order corrections, 
even though it may be a rather challenging task. 

One important question that we have left for future work is to study these effects 
from the point of view of the dual theory. Of course, this theory will be strongly cou- 
pled, and for general dimensions it is not clear what the precise dual theory even is. 
We can, however, construct phenomenological dual matrix models along the lines of 
[23\ [Ml [26l [27] and compare the various phase structures of the bulk theories that we 
have uncovered with those of the dual matrix models. It is also of interest to perform 
our analysis in the canonical ensemble where the charge, rather than the electric po- 
tential, is held fixed. In this case, the reference background will no longer be thermal 
anti-de Sitter space, but an extremal black hole. Since we have uncovered new black 
hole solutions, it is important to study their dynamical stability with respect to small 
perturbations and see whether the regions of negative specific heat and dynamical 
instability coincide [SH [23 [5S1 [SZ]- Finally, there has been a recent interest in 
higher derivative corrections to the shear viscosity bound and the Gauss-Bonnet-F^ 
black hole, and further higher derivative corrections thereof (particularly RF^ terms) 
may be of interest to this story [581 EH [60l [6T]. 
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A Equations of Motion 

Gauss-Bonnet Gravity 

The general action of gravity with matter is given by: 



/ = 



/ 




X^/^ (^R -A + aR^ + hR^^Rf"" + cR^upaR^""^"] + Imatter- (A.l) 



We are particularly interested in the Gauss-Bonnet combination. 



a = c. 



b = -4c. 



(A.2) 
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The equations of motion for the Gauss-Bonnet combination are given by: 

+ c {2RR^, - m^pRi - AR^p.aR"' + 2RlP'^R,xpa) ■ (A.3) 

B Useful Intermediate Results 

In order to make hfe simpler for one to verify the equations of motion with our ansatz, 
we give some of the intermediate steps of the calculation. 

ChristofFel Symbols 

The nonzero Christoffel symbols obtained from our ansatz are: 



yt 
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= V , 


■pr 


pr 
rr 


= A', 


n. 

IJ 


r* 


_^ 
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pi 



jk- 

Curvature Tensors 

The nonzero elements of the Riemann tensor: 

R,^^ = -e^-Rl,, = [v" - \'u' + v''') 
Rtitj = -e^^Rltj = r,y'e^^''-^^gij 

Rrirj = gij 

Rijki = f'^Rijki — r^e~'^^ {gikgji — gugjk) ■ 

The elements of the Ricci tensor are: 

R^^ = _ (j," _ AV + u'^) + ^^^^^ 

R^ = e2(-^) (u" - X'u' + v'^ + ^^^^) 

Rij = \k- [r [y' - A') + d - 2] } ~gij. 

Finally, the Ricci scalar is given by: 



(B.l) 



(B.2) 



(B.3) 
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+ 



+ (f^. (B.4) 



Gauss-Bonnet terms 

We list various relevant contractions for the various derivations. 
Quantities in the Action: 



2 [u" - AV + u'^)^ + 2 (d - 1) {u" - XV + u'^) ^'^ 



+ ^ {k - [r - A') +id-2)]y (B.5) 



R^.uXpR'"'^'' = 4e-^^ {i^" - AV + + 4 (d - 1) e 



-4A 



A'2 + I.'2 



2(d-l) (d-2) / _ 
^ ^ 'd-2 



(B.6) 



- AR.^R^'^ + R,.,,R^^'^ = - 1) I (^L^2)r4^^' 



+ e" 



-2A 



^ fc(i/^^- AV + i^^2) ^ ^ (rf-3)fc(i/^-A0 _ (d-4)fc 



+ (d - 2) e 



-4A 



4 0^^^ + '^^--3A^/) (d-3)(;/-A0 (r/ - 3) (d - 4) 



■• (B.7) 



R^pR'^P non-vanishing terms: 



l\ 2 



r\ 2 



(B.8) 
(B.9) 



e - 2e-'^^k {r {u' - A') + (d - 2)) 
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+e-^^ (r^ {u' - X'Y + 2{d-2)r {u' - A') + {d- 2)^)] g'K (B.IO) 



Ri^p'^'^ Rp^ non- vanishing terms: 



-2v-A\ 



+ 



i/^ (z/^ - X') _ {d - 2) u' ^ ^2xt!^ 

^3 ^3 



(B.ll) 



i^'-'^'-'^i?/,, = ((i - 1) e" 



-6A 



zy" - A'// + ;/2r z/' //'-AV + ;/ 



+ 



X' ju' - X') {d - 2) A^ 



(B.12) 



-4A 



(d-2)(i/-A') (d-2r 



-2A 



-k{r {u'-X')+2{d-2))~g'^ + -^g'^. 

(B.13) 



R^ K\pR^'^^'' non-vanishing terms: 



(B.14) 
(B.15) 
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Electromagnetic Intermediate Quantities 

We give some intermediate steps in the calculation of the matter energy tensor and 
the solution of the Maxwell equations: 



F^,Ff'^= -2f{ry 



(B.17) 



F^^F'''= -e 



pip FpcjF'^'^ Ft* 



e f[r) 

-2A ^ 



All other relevant terms are vanishing. 



-{u+\) 



fir) 



C Thermodynamic Derivatives 



(B.18) 

(B.19) 
(B.20) 



In this appendix we present the thermodynamic derivatives that allow us to compute 
the thermodynamic charge, entropy and mass used throughout the paper. Let us 
begin with the expressions in the rn, Q parametrization. For the grand canonical 
ensemble we used that: 



Q 



E 



1 
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{do) 
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( _d±\ _ ( d±\ \S^hJq _ 



(C.l) 



s = p 
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I, (C.2) 



(C.3) 
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The expressions for the specific heat and electrical permittivity are given by: 



= { ^ ) = ^ T^'F^- (C.5) 

^Th f9^\ _(d4L\ ^ 



and 



The above expressions can also be obtained in the rn, Th parametrization for which 
we have provided an thermodynamic analysis for the £ = case in appendix [Dl 

di 

Q = -Thj^, (C.6) 



S = -Th 



/, (C.7) 



E = ^-^ + Q(t>. (C.8) 

Furthermore, we provide the expressions for the specific heat and isothermal electric 
permittivity in the rn, Th parametrization: 



D The £ = case in th-, Th parametrization 

Even though we discussed throughout the text the practical usefulness of the rn, Q 
parametrization in analyzing the thermodynamics for non-zero £ due to the highly 
non-trivial dependence on Q in our expressions, this is not the case when £ = 0. It 
turns out that the analytic expressions we can obtain in this case are much simpler 
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in the rn, Th parametrization for e = and it is instructive to include the main 
results for the global and local thermodynamics. We note that the analysis in this 
parametrization for e = has been performed in [27^ I31j . This provides us with a 
non-trivial check of our interpretations of the e = analysis and further support for 
our interpretation of the e 7^ results. 

First we need to solve for ^ and Q in terms of rn and Th- We get: 



d-2 



d-A 



-AttTh 



{d-2Y 
(d- 



l)rH 



d-l 



+ 
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2kDrH 
J 



d-3- 



, (D.l) 



2r 



d-2 
H 



9^ 



{d - 1) hrn"^-^ 



{d - 4) Dh'^rn'^-^ 



{d-2) 



id-l)rH 



d-l 



+ 



2kDrH 



d-3 



d-2 

We find that the condition for a massless black hole can be written as: 



rpM=0 



{d - 2) K^Arjj + d{d-2)krjj + d{d- ck 
2TTd (d - 2) [r% + 2 {d-3) ckn'^rn] 



2,^2 



(D.2) 



(D.3) 



One important thing here is that we need to ensure that is always positive, so 
that the charge that enters our Lagrangian is real. In the rn, Q parametrization this 
was trivial; however, now it is not. The payback in our parametrization is that the 
extremality bound has become trivial. We find that Th needs to satisfy: 



Th < 



{d - 2) K^Arj^ + {d-l){d- 2) kr'jj + {d - I) {d - 3) {d - 4) ck'^n 



4tt {d -l){d- 2) [r% + 2 {d-3) ckn'^rn] 
Finally, we substitute and Q in the free energy to find: 



(D.4) 
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1) {d - 3) k'cr'lj 
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, (D.5) 
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from which we can obtain the following critical temperature: 

A 4 (d-l)(d-3)fc^c 

4^ ^-fcc(d-l)rH 

Finally, we provide the expressions for the specific heat and isothermal electrical 
permittivity of our system. The thermodynamic derivative expressions for these 
quantities are given in appendix O The specific heat at constant electric potential is 
given by: 



_ 87r2 {d - ly r'jf^TH {rjj + 2c{d-3)k 



k2 



2 (d - 1) ttthTh {r% -2c{d- 3) kK^) + (^c^J^^^^^§^k^ + Arj^ 

(D.7) 

The electrical permittivity at constant temperature is given by: 



g 



Sfc ~ 2 (d - 2) 

2^rH ((2d - 3) 4 + (2d - 5) fc) - T (A; ((d - 2) r2^ + iitgl^^^ _ ^2Ar-^ 

(D.8) 

The above expressions give rise to critical temperatures for the thermal and electrical 
local stability of our thermodynamic system. The critical temperatures indicate a 
change in sign of the specific heat and electrical permittivity and thus a change in 
stability. The critical temperature for thermal stability is given by: 

27r (d - 1) (rl, - 2 (d - 3) ckn^) m ' ^ ' 

The critical temperature for electrical stability is given by: 

27r (d - 1) {rjj - 2 (d - 3) ckK^) m ' ^ ' 

(d-l)fc 2 I c{d^3f{d-4) j2 A„4 
r(2) = d-2 ^ ^^^H 

27rrH (^r2^ + 2 (d - 3) (2d - 5) cfc) ■ ^ ' ' 

Finally, we give the diagrams corresponding to those discussed in section [5l 
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Fig. 25: Thermal structure for /c = 2, d = 4 and (a) c < g^l|^| , (b) c > 
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Fig. 26: Thermal structure for d > 4, k = — (d — 2) . The two regions of c displayed 
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Fig. 27: Thermal structure for d > 4, = — (d — 2). The two regions of c displayed 
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E Explicit regions 



E.l c-(j? Plane 

Before proceeding to calculate the form of the curves separating regions with qualita- 
tively different behavior in the c-cjP' plane, it is useful to write down the temperature 
as function of the horizon radius and the electrostatic potential, thus making it simple 
to have the form of a constant (j) slice in the rj^, Th plane. Using equations (|4.1U|) 
and ()4.16p we easily find that: 



A; > 

First of all, let's study the region c > 0. We can see in figures [2]^ a), (b) andlU that 
a constant (p curve does not cross the F-curve if (j) > (j)o, where (po is the potential 
for the curve that incudes the extremal black hole with vanishing free energy. So 
requiring that Tc {rE,F=o) = we can find the radius of this black hole: 



4 _ did-l)id-3)ck^ 



So solving T {rE,F=o-, 0o) = we get: 




The above curve obviously crosses the cj) axis at (pQ = Notice that for d > 4 

a constant (p curve with <p > (pQ, has two disconnected branches, one where the free 
energy is always positive, and one where the free energy is always negative, while for 
d = 4 there is only one branch where the free energy is always negative. 



Let's now study the temperature behavior as function of the potential. We can 
see in figures [JJa), (b) and HI that all curves with <p > (pQ, where (po is some critical 
potential, contain extremal black holes of finite size. Actually, for d = 4 these curves 
contain exactly one extremal black hole, while for d > 4 they contain two. These two 
can be either both globally unstable, or one stable (the larger one) and one unstable. 
In order to specify this critical (po, we need to find the radii of the extremal black 
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holes solving T {vh, (p) = 0. This can be written as: 

_ + ( - 2 - 2f A rl + (■'-l)(j-3)(^-4) ^,. ^ „^ 

\ K'^ J a — I 

This is just a quadratic equation for r"^. We can always check if there are real roots 
using the discriminant. Obviously (j)Q is the potential that sets the discriminant equal 
to zero. That leads us to: 



a2 



ii-l)k (.-l)(.-3)(d-4) 

° 2 ((i- 2)2^2 V d{d-2f ^ ' 

Notice that if there are real roots, then they are both positive, as we noticed in figure 
m Notice also that for d = 4, the above curve reduces to a constant 4> curve. Finally, 
it crosses the axis at the same point as the previous curve. 

Let's now consider the region c < 0. Observing figures [3l [5]^a) and (b), [6||a) and 
(b) and [71 we see that a constant ^ curve enters the negative entropy region, if 
(f) < (pQ where (po is the potential for the curve that contains the extremal black hole 
with zero entropy. From equation (j4.14p . we know that for ck < all black holes 
with vanishing entropy have horizon radius: 

rs=o = \/-2{d-l)ckK^. (E.6) 
So requiring Th {rs=Q-, (po) = we find that cp^ equals: 

" A(d-2fK^ (d2 + d - 8 + 4 (d - 1) ((i - 2) CK^K) . (E.7) 



Now in order to check whether a constant (p curve crosses the local stability curve, 
first we check the behavior of the temperature as r — > r^i^, where 



r^nin = V-2(d-3)cfcK2. (E.8) 
We solve lim T (r, 0o) {''' — fmin) = 0, to find: 

T mill 

= ^o^3 2 {d{d-l)+iid- 2) (d - 3) ck'^A) . (E.9) 
4 (a — 2) 
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The above curve crosses the (j) 'ixis at (/)q = ^^~^y^~2' which is lower than the point 
where the curves lying in the c > region cross the (j) axis, except for d = 4 when 
they coincide. For potentials smaller than the above there are two black holes for 
every temperature, and for larger one or three black holes, except for d = A when 
there is always one. 

Notice that the behavior of a constant (p curve as r ^ rmin determines whether 
this curve crosses the local stability curve an even or odd number of times. In[2]^a) 
and [5^a) and (b), we see that there are some cases when the local stability behav- 
ior changes without change in the behavior of the temperature for t — > ^min* 

This 

happens for positive c at d = 4 and for negative c at d > 4. In both cases the local 
stability curve has a minimum temperature, and there are three black holes for the 
same temperature (for some temperatures) \i (p < (pQ, where (j)Q is the potential for the 
curve that contains the aforementioned minimum. So in order to specify the curve 
that separates the area with one black hole from the area of three black holes we need 
to solve 

dTc, (r) 



0, (E.IO) 



dr 

T (rininT,C^=0, </'o) = (?^mm T,C^=o) • (E-H) 

For d > 4 this gives a quite complicated expression. Some important characteristics 
of this curve is that it crosses the curve separating the odd-even number of black holes 

areas at the point = ~ 4(^d-i){d-2)K'''A '' '^o ~ ' ^'^^ ^ smaller than the 

above the local stability curve is one-to-one in rHjTn and there is no constant (p 
curve containing three black holes with the same temperature. Finally, it crosses the 
(j) axis at the same point as the curves at the c > region. 

Finally, at d = 4, c > 0, it is easy to solve the relevant equations to get 

^mmT,C^=0 = QckK^ , (E.12) 



,1 3^ 



^ [l + 6ckK^) . (E.13) 

This curve crosses the c axis at c = — . There are no constant (p curves containing 
three black holes at the same temperature for c greater than the above. 

Let's check a detail about the d > 4, c < case. It is obvious that for some potentials 



62 



the metastable black hole (the smallest one) has manifestly negative entropy. This 
happens for < 00, where 0o is the potential for the curve crossing through the black 
hole with infinite specific heat and zero entropy. Solving Tc^ {rs=o) = T {rs=o, 0o) 
we get 



k 



4(^-2)^^2 



{(f - + 4d - 8 + 4 (d - 1) (d - 2) (d - 4) ck^A) . (E.14) 



The above curve is tangent at the curve separating the one black hole region from 
the three black holes region at c = - 4(fe^^^2®^5)7?A • 

Finally, let's be a little careful about the Hawking-Page phase transition for c < 0. 
It is possible that a constant (f) curve crosses the global stability curve; however, 
this happens at an area of locally unstable black holes or negative entropy black 
holes. The first case occurs if c > ~ 4(d-t)td^l)K'^A 4* > 'Poi where 0o is the 
potential for the curve crossing the mutual point of the global stability and local 
stability curve. We can specify this point by solving = 0, and then 

»'='"minT,F=0 

T {rYnmT,F=o, (po) = ('^minT,F=o)) to get a very complicated expression. 

The second case occurs when c < ~ 4{d-'i){d-2)K-^A 4* > 4'0 where (po is the 
potential for the curve crossing the mutual point of the global stability and zero 
entropy curve. We need to solve Tp=o {rs=o) = T {rs=o, (po) to get 

= 2 {d{d+l)+^{d- 1) {d - 2) ck'A) . (E.15) 

4a [a — 2) 

The two curves above meet at = ~ 4(d-i)^(7-2)w^A ' 4o = 2{d-2fK^ ) ' ^ 
they are also parallel, thus forming a continuous and smooth curve. For c equal to 
the above the entropy bound occurs at exactly the minimum of the global stability 
curve. 

One final detail. If one extends the curves above, in the complementary regions 
of c we used, there is an area created, containing curves with three Hawking-Page 
phase transitions. However, always only one of them happens in a region with both 
positive entropy and positive specific heat. Specifically in the relevant region for 

c > ~ 4(d-i j^rf-2)K"tA ' black holes smaller in radius, have negative specific 

heat, and the smallest has negative entropy, while in the case c < ~ 4(d-i)^d^2)K-^A ' 
the two smaller ones have negative entropy and the smallest has negative specific heat. 
In this region the curve separating odd-even number of black holes appears, and once 
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one enters in the region of two black holes, there are actually only two positions of 
Hawking-Page phase transition, again only one being valid, as the other is in a region 
with negative entropy. This happens for c as small as c = ^fi_2)^d~l^2d-l)K-^A ' '^'^^^e 
our curve ends on the odd-even separating curve. 

k<0 

Observing figuresEI^a), (b) and[9]|^a), we see that a constant <p curve enters the negative 
mass region (p < (j)Q, where (po is the potential for the curve that contains the 
extremal massless black hole. We can find the horizon radius of the extremal massless 
black hole using equation (ID.Sh : 



So requiring that Th {rE,M=o, (po) = 0) we find that cpQ equals: 



»S = ^ 1 + iSi^l. (E,7) 



'd{d-2)' 

This curve crosses the c axis at c = ~ i{d-i){d-l)K'^A • ^ larger than the above 
there are no negative mass black holes. 

Observing figures [8l^b) , [9|^a) and (b) we see that a constant 4> curve enters the neg- 
ative entropy region, if </> < where 0o is the potential for the curve that contains 
the extremal black hole with zero entropy. From equation (j4.14p . we know that for 
cA; < all black holes with vanishing entropy have horizon radius: 



rs=o = \/-2{d-\)ckK^. (E.18) 
So requiring Tu (rs=^^ 0o) = we find that </)o equals: 

" ^{d-2fK^ (d^ + d - 8 + 4 (d - 1) (d - 2) CK^A) . (E.19) 

The above curve crosses the c axis at c = ~ A[d-\){d-2)K'^K ■ This is the minimum c for 
which there are negative entropy black holes. 
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Now trying to combine the above, we check whether a constant (j) curve enters a 
region of negative mass black holes with positive entropy. If we observe figure WL^) 
this happens if < 0O) where 0o is the potential for the curve that contains the 
massless black hole with zero entropy. Thus solving Th {rs=o^4>o) = T^^^ {rs=o), 
we get: 

(d(d+l) + 4((i-l)(d-2)cK^A). (E.20) 



/,2 ^ ^ 

° 4d(d-2)2-2 



This curve crosses the massless curve and the entropy curve, at the same point 

= - 4id-}nd-l)n^A ^^0 = -JTTjih^)- This value of c is the minimum for 
which there are positive mass black holes with negative entropy. This curve also 
crosses the c axis at c = — 4^^_^j^^2)K^A • This is the maximum c for which negative 
mass black holes with positive entropy exist. 

Finally, in order to know if a constant (j) curve crosses the local stability curve we 
need just to specify the behavior of the temperature as horizon radius approaches the 
minimum radius 



Train = \/-2 (d - 3) cfe^^ (E.21) 

as we can see in figure El^b). If it the curve goes to infinity, that means that it had 
a minimum at some horizon radius. There it has crossed the local stability curve. If 
is goes to minus infinity then from the figure we can see that it is one-to-one, and 
obviously it has not crossed the local stability curve. So the critical (pQ is the one for 
which lim T (r, (po) (r — r^i^) = 0. Solving this we get 

^ {d{d-l)+4{d-2){d-3)cK'^A) . (E.22) 



° 4 (d - 2f 

This curve crosses the c axis at c = — ^^^_^j^^g.^^4^ . This is the minimum c for which 
there are locally unstable black holes. It coincides with the minimum c for which 
there are no negative mass black holes. 

E.2 e-(j)-^ Plane 
A; > 

For e < the e-bound excludes the region of the black holes where the small x 
expansion is good. That allows us to specify the regions of qualitatively different 
thermal behavior approximately, using the small x expansion. 
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First let's specify which constant (j) curves contain an extremal black hole. We can 
see from figure HD^a) , that extremal black holes are contained in constant (j) curves 
only \i (j) > 4>E, where 0s is the potential of the curve containing the extremal black 
hole that saturates the e-bound. Using the zeroth order forms for the potential ()4.10p 
and the temperature (|4.16p . we acquire the result: 



{d-l)k 



{d-2)K\l 2 (-108^4 Ae + 1) ' 
while if one uses the first order corrected ones ()7.5p and ()7.6p . one finds: 



(E.23) 



85(i-116 / {d-l)k 



^ 3{d- 2) (3(i - 4) K y 6 (-2916^4 Ae + 29) ' 
Notice that both results are of the form: 

A 



(E.24) 



(E.25) 



e — El 



and actually are very close to the actual numbers we calculated numerically, taking 
into account the full dependence on e. 

Similarly in order to calculate which potentials cross the F-curve, thus allowing 
for Hawking-Page phase transitions, we need to specify the potential containing the 
black hole with zero free energy that saturates the e-bound, (j)HP- Then a constant </> 
curve crosses the F-curve provided that cp < (j)HP- Once again using the zeroth order 
result for the potential and the free energy we get: 



1 / d{d-l)k 

VHP ~ 



(d - 2) K y 2 (108 {d - 2) g^Ae + d) ' 
while at first order: 

85d-116 



(E.26) 



VHP 



d{d-l)k 



3 (d - 2) (3d - 4) K ^ 6 (^2916 {d - 2) g^Ae + ^ 



(j(83d-112) 
4 



(E.27) 



We again observe that both results are of the form: 

V^- (e - eo) 



(E.28) 
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but now the value of £q depends on the number of dimensions. It looks like Eq will 
equal —ef at the limit of infinite number of dimensions. 

We can specify similarly the curve separating regions with one and two black holes 
to find at zeroth order: 



1 {d-l)k , , 

'^^^= (d^V2WA. + l)' ^^-^^^ 

Notice that this curve is of the same form as the curve specifying the existence of 
extremal black holes, where we have substituted e with —e. Thus the positions of 
the poles of the two curves are opposite, and so are the inclinations of the two curves 
in the e-cj)"'^ plane. This agrees very well with our numerical results, where we have 
taken into account the full e dependence. 

We further notice that the inclination of the curve separating Hawking-Page and 
non-Hawking-Page constant (f) curves in the e — ^ plane is a fraction of that of 
the curve separating one and two black hole solutions, the later being equal to the 
opposite of the curve separating constant (p curves containing extremal black holes. 
This is consistent with our numerical results. 

We can specify specify the values of ef and in figure [T71 is the value of e 
that the extremality curve crosses the axis, thus it is the minimum e, that extremal 
black holes exist. Using the first order result we acquired above we find: 

(E.30) 



1 2916C/4A' 

Finally, is the value of e where the Hawking-Page curve crosses the extremality 
curve and ends. Again using the first order results we got above we find: 

= 2916 (3d - 4) 5^ A- ^^-^^^ 

Unfortunately for e > the critical constant (p curves cross the F-curve and C^- 
curve in regions where the small x expansion is not good, thus making the numerical 
calculation of the different regions necessary. 

k<0 

As we can see in figure [T6l a) and (b) for e < 0, the constant cp curves are again 
separated in two classes by the curve containing the extremal black hole that saturates 
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the e-bound. Repeating the same exact procedure as in the A; > we find the same 
exact result given by equations (|E.29p and (|E.24p . Notice though that now the curves 
containing extremal black holes are those with (p < (j)E- 

Finally, we can calculate which constant curves contain massless black holes. 
They are those with (p < (pM, where (pM is the potential of the curve containing 
the massless black hole that saturates the e-bound. Using the expressions for the 
thermodynamic energy (j7.1ip we get at zero order: 



d{d-l)k 



2 {d - 2) {108 {d - 2) g^Ae + d) 



(E.32) 



and at first order in e: 



85d - 116 

(PM 



d{d-l)k 



3 - 2) (3d - 4) K ^ 6 (^2916 {d - 2) g^Ae + ^^^f^) ' 



(E.33) 



Notice that strangely these expression are identical with those for (pHP for the k > 
case, having substituted k with —k. The inclination of this curve in the e — ^ plane 
is a fraction of that of the curve separating regions with extremal and without 
extremal black holes. 

Let's specify the values and in figure [HI is the value of e that the 
extremality curve crosses the axis. As in the k > case we find: 

99 

r+ (E.34) 



29I654A 



ej~, is the value where the curve separating constant (p curves containing and non- 
containing massless black holes meets the extremality curve and ends. Using the first 
order results we acquired above we find: 

d(85d-lU) 
- ^ ' (E.35) 



2916 (3(i- 4) g'^A' 



Once again we cannot acquire the critical curves using the small x expansion in 
the e > case. 
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